Optical Properties and the Electro-optic and Photoelastic 
Effects in Crystals Expressed in Tensor Form 

By W. P. MASON 

I. Introduction 

THE electro-optic and photoelastic effects in crystals were first investi- 
gated by Pockels, 1 who developed a phenomenological theory for these 
effects and measured the constants for a number of crystals. Since then not 
much work has been done on the subject till the very large electro-optic 
effects were discovered in two tetragonal crystals ammonium dihydrogen 
phosphate (ADP) and potassium dihydrogen phosphate (KDP). With these 
crystals light modulators can be obtained which work on voltages of 2000 
volts or less. Their use has been suggested 2 in such equipment as light valves 
for sound on film recording and in television systems. Furthermore, since 
the electro-optic effect depends on a change in the dielectric constant with 
voltage, and the dielectric constant is known to follow the field up to 10 10 
cycles, it is obvious that this effect can be used to produce very short light 
pulses which may be of interest for physical investigations and for strobo- 
scope instruments of very high resolution. Hence these crystals renew an 
interest in the electro-optic effect. 

In looking over the literature on the electro-optic effect and photoelastic 
effect in crystals, there do not seem to be any derivations that give them 
in terms of thermodynamic potentials, which allow one to investigate the 
condition under which equalities occur between the various electro-optic 
and photoelastic constants. Hence it is the purpose of this paper to give such 
a derivation. Another object is to give a derivation of Maxwell's equations 
in tensor form, and to apply them to the derivation of the Fresnel ellipsoid. 
The first sections deal with the optics of crystals, and derive the Fresnel 
ellipsoid from Maxwell's equations. Other sections give a derivation of the 
two effects, discuss methods for measuring them by determining the bi- 
refrigence in various directions and give the constants for the two effects in 
terms of crystal symmetries. The final section discusses the application of 
the photoelastic effect for measuring strains in isotropic media. 

1 F. Pockels, Lehrbuck Dcr Kristalloptic, 13. Teubner, Leipzig, 1906. 

2 See Patent 2,467,325 issued to the writer; "Light Modulation by P type Crystals," 
George D. Gotschall, Jour. Soc. Motion Picture Engineers, Julv, 1948, pp. 13-20: B. H. 
Hillings, Jour. Opt. Soc. Am., 39, 797, 802 (1949). 
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II. Solution of Maxwell's Equations In Tensor Form 
In tensor notation, Maxwell's equations for a nonmagnetic medium with 
no free charges take the form 

1 dDi _ dHj 1 djh = _ dEk . dDj = dHj = (1) 

V~dT~ Uik dx k ' V dt ' lkl dXi ' dxt ' dx,- 

where Di is the electric displacement, Hy the magnetic field, E k the electric 
field, V the velocity of light in vacuo and e ijk a tensor equal to zero when 
i = j or k orj = k, but equal to 1 or — 1 when all three numbers are different. 
If the numbers are in rotation, i.e. 1, 2, 3; 2, 3, 1; 3, 1, 2 the value is +1 
while, if they are out of rotation, the value is —1. 

We assume the electric vector to be representable by a plane wave whose 
planes of equal phase are taken normal to the unit vector m . Then 

E h = E J uit - XiniM (2) 

where E 0k are constants representing the maximum values of the field along 
the three rectangular coordinates and j = \/^\. Substituting (2) in the 
second of equations (1), noting that E 0k are not functions of the space co- 
ordinates, we have 

1 dHj J03 r p 1 Jmlt-mmM (7) 

v~af = H ^ jhih0kUii * ■ ^ } 

Integrating with respect to the time 

II j = - \(jkiEn k ni\e — tia.e . v*) 

v 

Hence, 

Ha; = - UjkiEo h nj\ (5) 

' v 

and therefore the magnetic vector is normal to the plane determined by 
Eo k and w, . 

Next, using the first of equations (1), 

= — — [fin, H o; nic\e 

V 

Integrating with respect to time, 

Di = -^[e ijk H 0i n k ]e M '- Ik " kM . (7) 

v 
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Inserting the value of H 0i from (5), this equation takes the form 

i= r Uuk \€jki •Eo k niJtikle 

1r 



V 2 



and, in general, 

Di = —^■Uijk(ejkiEi c Ui)nic]. (9) 

v 

Expanding the inner parenthesis, we have the components 

(Eoiii — Em)i] (E 3 iii — Ei>h) 2 ; (Eiii 2 — Eoni) 3 . (10) 

Then 

€ t[(£ 2 »3 - E-jh); (Emii - Ei)i 3 ); {E x n 2 - E 2 ih)\n k gives 

V 2 
A = 7 {{E-iHi — Einz)n 3 — (Ei>h — E 2 m)n 2 ] 

v- 

= [(E 3 n 3 + E 2 Mo + Ei m)ni — Ei(nl + n\ + »*)] 



V 
D> = — — r [(£i«-> — Eom)ni — (/v-Mt — £3»2)»s] 
a 2 



(ID 



= [(£3^3 + E 2 n> + Eim)n 2 — E 2 (m + m 2 + W3)] 

Z) 3 = — — -[(£ 2 «j — E 3 n 2 )n 2 — {Eon — Eim)ni\ 
v i 

= [(£ 3 "3 + £2 "2 + -Ei "i)«3 — £:t(«i + n 2 2 + »!)]. 
Now, since n\ + »! + »* = 1 because n is a unit vector, we have 

Di = — [JS< - (Ejnjhi] or -^ Z>, - £, - (£,»,)»< = 0. (12) 
v- V 

This equation states that Di , £, and m are in the same plane, Hj being 
normal to the plane as shown by Fig. 1. The energy flow vector 

Si = ] —eij k EjH k (13) 

• 47T 

also lies in the plane since it is perpendicular to E and H. It is at the same 
angle with n that E is with D. The velocity of energy flow is v/cos 6. The 
energy velocity is called the ray velocity and the energy path the ray path. 
Next, from the relation for a material medium, that 

Di = KijEj or conversely Ej = fijiDi (14) 
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where Ka are the dielectric constants measured at optical frequencies and 
(jji are the impermeability constants determined from the relations 

fe = A y, '/A A ' (15) 

where 

Ku K\i K\3 

A = K12 Km K23 

K\s Ka K33 

and A" the determinant obtained by suppressing the/ h row and i tl column, 
we can eliminate E t from equation (12) and obtain 
2 
^zDi - 0n A + /3 12 A + /3 13 A - {EjnDm 

2 
^5 2 = P12 Di + P22D2 + 023 A - {Ejn,)n2 (16) 

2 
^ A = fiuDi + 023 D 2 + /J33A - (EjnM. 

This can be put in the form 

(£/»>i - ADSii ~ &/V 2 ] + J812A + ftsA 
(£^)«9 = &2A + 0^22 - */F*)A + ft«D 8 (17) 

(£^y)« 3 = /3uA + /3 23 A + 03 33 - dV^A- 
Solving for A , A and A 

A = [05« - ttyFWss - */P) - /sSjlfiiffiK 
A - [03 u - ^/FOfos - */P) - /3i 3 ][JSy»/]«a (18) 

A = [03ii - */V*)</8n - ^ 2 /n - PiaPi^K 
Now, since D and « are at right angles, 

Awi + A«a + A«s - 0. (19) 

Hence, 
= [(0 22 - * 2 /F)63 33 - *»/P) - /&]»? 

+ [(fin - »»/F)(j8„ - ^/F 2 ) - |&]«f (20) 



+ [Cfci - ^ 2 /P)03 2 2 - z^/F) - /3 2 2 ]» 3 . 
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Fig. 1 — Position of electric, magnetic and normal vectors for an electromagnetic plane 

wave in a crystal. 




*-Y 



Fig. 2 — Rotated axes and angles for relating them (o unrotated axes. 



By choosing the original v, y, z axes so that ffu = 0i 3 = /3 23 = and using 
the values /3u = /3i , /$> 2 = #> , /3 33 = 3 this gives the equation 



2— 9 + — ^ + -^-, = 0. 



v ~ v „ v 

ft " Fi A - F ft - yi 



(21; 



For transmission along the .Y axis «i = 1, «j = >/ 3 = and the two velocities 
are given by 



& = p 2 v°- = b°; 1? = /3 3 F 2 = c 2 . 



(22) 
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Similarly the third velocity v- = faV 2 = a 2 can also be used and equation 
(21) reduces to 

n * + w » + w ' = 0. (23) 

a* - v* ^ b* - & c 2 - v 2 

This is a quadratic equation for the velocities v in terms of the principal 
velocities a, b and c which are usually taken so that a> b > c. 
Solving for the velocities, we obtain the quadratic equation 

*■' - i?[n\(b 2 + c 2 ) + n\{a 2 + c 2 ) + 4(a 2 + & 2 )] 

, •> 2 * ' 

+ nlb 2 c 2 + thfP& + n 3 a 2 b 2 = 0. 

Letting L = n\{b 2 - c 2 ), M = n\(c 2 - a 2 ), N = ih(a 2 - b 2 ) the solutions 
for the velocities become 

2v 2 = n\(b 2 + c) + n\(c + a) + nl(a + 6 2 ) 

± y/l? + M 2 + N 2 - 2LM - 2LN - 2MN 

This equation can be put into a simpler form if we change to the coordinate 
system shown by Fig. 2. Here the rotated system is related to the original 
system by three angles 9, <p, \f/. 6 is the angle between the Z axis and the 
Z axis, ip is the angle the plane containing Z and Z' makes with the X axis 
while \f/ represents a rotation of the primed coordinate systems about the 
Z' axis. The direction cosines for the primed system with respect to the 
normal system are designated by the matrix 

X Y Z 
X' t x mi m (26) 

U m 3 n 3 



y 

Z' 



where, in terms of 6, <p and \f/, these direction cosines are, 

/i = cos 6 cos <p cos if/ — sin <p sin \p, 

nil = cos 6 sin <p cos \p + cos ip sin \p, rt\ = — sin cos \p 

f 2 = — cos cos ip sin i/< — sin <p cos ^, 

tn 2 = cos <p cos ^ — sin ip sin ^ cos 0, n 2 = sin 5 sin \J/ 

U = cos <p sin 0, m-i = sin ip sin 5, » 3 = cos 0. (27) 

If we take Z' as the direction of the wave normal, then in equation (25) 

»1 = (3 , >h = »3 , "3 = »3 
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and the equation for the velocities becomes 

2v ■■= a 2 (sin 2 (p sin 6 -f cos"0) + b~(cos~ ip sin" 9 -\- cos"0) + c sin 2 9 
2 — 6 8 ) 2 (cos 2 9 cos- <p -f- sin 2 ip) 2 + 2 (a 2 — b 2 ){c 2 — b 2 ) 
sin 2 0(cos 2 9 cos 2 <p — sin 2 ip) + (c 2 — & 2 ) 2 sin 4 q 



J* 



(28) 



A very elegant construction for the wave-velocities and the directions of 
vibration is the Fresnel index ellipsoid. Consider the ellipsoid 

fl 2 T 2 + b y„ + C 2 2 2 = i (2Q) 

Then Fresnel 3 showed that, for any diametral plane perpendicular to the 
wave normal, the two principal axes of the ellipse were the directions of the 
two permitted vibrations, while the wave velocities were the reciprocals of 
the principal semi-axes. 

We wish to show now that the maximum and minimum values of the im- 
permeability constants in a plane perpendicular to the direction of the 
wave normal determine the directions of vibration and the values of the two 
velocities. To show this we make use of the fact that /3,-y is a second rank 
tensor and transforms according to the tensor transformation formula 

fi'„ = P £' ft, (30) 

dxk dxt 

where the partial derivatives are the direction cosines 



dxi 

dxi 


dx[ 

d.t 2 


dxi 

— - = Mi 

dx 3 


dXi . 
dx\ 


dxj 

= Wo, 

dx 2 


dxi 

— = n- 

dx 3 


dx 3 

t '3, 

dxi 


dxs 

a,, = "•■ 


dXz 

— = n 3 
dx 3 



Expanding equation (30) the six transformation equations become 
ft'i - tffri + 2timi0u + 2tw$n + m\fa + 2*» 1 » J /3 23 + n%* 
fi'n = &60u + {dm + mfafa + {tm + *i4)0u + miflfe/fe 

+ (»h>h + Mim; 2 )/3 2 3 + nahPM 

/3i3 = (lftfu + (Aw's + »hh)0i» + (A»a + >hQ@n + mmsfoi 

+ (»iw : , -+- nii)h)l3^ + »i»a033 (31) 

3 Sec for example "Photoelasticity," Coker and Filon, Cambridge University Press, 
pages 17 and 18. 
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ft' 2 = 4fti + 24w 2 /3i2 + 24«2/3is + «4fe + 2w 2 w 2 ft 3 + w'fts 

023 = 66011 + (4»» 3 + «teWl2 + (4»3 + W2W1S + WfeWjftj 

+ (w 2 w 3 + «2fw 3 )ft 3 + «a»a/3sa 
/3 33 = $fci + 2/sf»a|8i2 + 2^i3/3 13 + wife + 2w 3 «3ft 3 + wsfts- 

Now, if the axes refer to the axes of a Fresnel ellipsoid, fe = /3i 3 = ft 3 = 
and one of the impermeability constants for any direction, say ft 3 , can be 
expressed in the form 

033 = 4ft + mjft + »2|3s (32) 

If r, which lies along Z' of Fig. 2, is the radius vector of the Fresnel ellipsoid, 
then the direction cosines 1% , w 3 and « 3 are 

x y z 

U = -, m 3 = -, «s = -• 
r r r 

From equation (24) ft = a 2 /F 2 , ft = 6 2 /F 2 , 0s = c 2 /V 2 and equation (32) 
becomes 

2 T r2„' 2 2 . t2 2 , 22 1 

fFft3=fl^ + by + cz =1. 

Hence the square of the radius vector of the Fresnel ellipsoid is 1/V fiat 
and the radius vector of the impermeability ellipsoid agrees with that of the 
Fresnel ellipsoid. Hence, the directions of vibration can be determined from 
the principal axes of the impermeability ellipsoid for any diametral plane. 

When light transmission occurs along Z', the direction for maximum and 
minimum impermeability can be obtained by evaluating fti and deter- 
mining the angle \J/ for which it has an extreme value. Inserting the direction 
cosines C\ , nil and »i from equation (27), we find 

/ „ T 2 „ 2 2 , sin 2m sin 2\f/ cos 6 , . 2 -2,1 
jSn = ft cos cos <p cos xp Z — _I + sin <p sin xf/ 

[■> „ . 2 2 , , sin 2<z> sin 2^ cos . 2 -2,1 (33) 
cos sin p cos if/ -\ ? h cos p sin \H 

+ /3 3 sin 2 cos i/\ 

Differentiating with respect to xf/ and setting the resultant derivative equal 
to zero, the value of xf/ that will satisfy the equation is given by 

(182 - ft) sin 2<p cos d 



tan 2\p = 



(ft ~ ft) (cos 2 cos 2 <p - sin 2 <p) + (ft - ft) sin 2 

. 2 2x (34) 

(b — a) sin 2<p cos 

(a 2 - b 2 ) (cos 2 cos 2 <p - sin 2 ^) + (c 2 - 6 2 ) sin 2 0' 
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For a given value on the right-hand side there are two values of \p, 90° apart, 
that will satisfy the equation and hence we have two directions of vibration 
at right angles to each other. Inserting (34) in (33) the values of fti and 
/3n for these two directions are 

20ii = /3i(sin 2 <p sin 2 0-f- cos 2 6) + ft (cos 2 <p sin 2 -f- cos 2 0) + 03 sin 2 



^ 



., - ft) 2 (cos 2 6 cos 2 <p + sin 2 vY- + 2(ft - ft) (ft - ft) 
* • sin 2 (cos 2 cos 2 ^ - sin 2 <p) + (ft - ft) 2 sin 4 0. 



Since ft corresponds to a 2 , etc., this equation agrees with the two velocities 
given in equation (28) and shows that the directions of vibration correspond 
with the maximum and minimum values of fti . 

It can also be shown that the two directions of electric displacement co- 
incide with the two values of \p given by equation (34). Transforming the 
electrical displacements to the X', V, Z' set of axes we have 

D [ = pDt + p D, +pD z = tiD t + m,D 2 + 1h D 3 
dxi ax-2 ox 3 

D[ = pD 1 + d ^D 2 + pD 3 = t 2 D, + m 2 D 2 + n 2 D 3 (35) 
dxi dx 2 Ox? 

D' 3 = pD 1 + pD 2 + pD 3 = l 3 D x + m 3 D 2 + n 3 D 3 . 
dxi ox 2 0X3 

Hence, inserting the values of D h D 2 , D 3 from equation (18), we find 

D[ = Wft - jSi'Ofo - /3ii) + w lW3 (/3i - /3ii)03 3 - Pu) 

+ nwztfi - /3n)(ft - ft'i) 

D 2 = Mfc - /3ii)(03 - ft'i) + »hm(Pi - 0u)(0, - ft'i) 

+ thmfa - ft'i)(ft - ft'i) 

D 3 = fl(j3 2 - ft',) (ft, - ft',) + »rM - ft\)(ft - jslo 

+ «iCSi - ft'i)j3 2 - ft'o. 

From equation (20) with ft* = ft:j = ft:s = 0, it is evident that the D 3 com- 
ponent vanishes and hence the two values of electric displacement lie in a 
plane perpendicular to Z'. By inserting the values of fti and the value of 
yp found from equation (34) we find that D» = and hence the electric dis- 
placement lies along the directions of the greatest value of fti . Similarly, 
from the second value of fti , D x vanishes and hence the second wave is per- 
pendicular to the first and in the direction of the smallest value of fti . 



(36) 
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III. Location of Optic Axes in a Crystal 

When the expression in the radical of equation (28) vanishes the two 
velocities are equal and an optic axis exists. Since the expression inside the 
radical can be written 

[(a 2 - 6 2 )(cos 2 0cos 2 v> + sinV) - (b 2 - c 2 )sin 2 0] 2 

- 4(a 2 - b 2 ){c 2 - b 2 ) sin 2 0sin 2 <? = 

then, since the square is always positive and since (a 2 — b 2 ) > and 
(b 2 — c 2 ) > 0, the equation can vanish only if ^ = 0. But <p = indicates 
that the two optic axes always lie in a plane perpendicular to the inter- 
mediate velocitv b. With <p = then the square vanishes when 



tan ! »=5k-g - 'an«= ±4 /f=|. (38) 

If (ar — b 2 ) < (b 2 — c 2 ) the value of the tan is less than unity and the 
crystal is called a positive crystal. For this case the two axes approach more 
closely the Z axis having the velocity c than they do the X axis. If 
(a 2 — b 2 ) > (b 2 — c 2 ) the crystal is negative. 

If a = b or b = c the crystal has a single optic axis and is respectively a 
positive or negative uniaxial crystal. For the first case the two velocities 
are given by 

i>i = a - b, v 2 = \/a 2 cos 2 d + c 2 sin 2 6. (39) 

The first velocity is that of the ordinary ray while that of the second is that 
of the extraordinary ray. Since a > c, the ordinary ray will have a velocity 
greater than the extraordinary ray except along the optic axis where they 
are equal. Since c < a, the maximum axis for any ellipse, formed by inter- 
secting the Fresnel ellipsoid at an angle to the optic axis, will lie in the plane 
formed by the normal and the c axis and hence the direction of polarization 
of the extraordinary ray will lie in the c, n plane. The polarization of the 
ordinary ray will be perpendicular to this plane. 

If b = c the a axis is the optic axis and the velocities of the two rays are 
again 

Vi = c and v\ = a 2 (l-sin 2 0cos 2 <p) -f- c 2 (sin 2 0cos 2 <p) (40) 

Hence, when 0= 90°, <p = 0°, the two velocities are equal and a is the optic 
axis. In this case the velocity of the extraordinary ray is greater than that 
of the ordinary ray except along the a axis, and the crystal is a negative 
uniaxial crystal. The polarization of the extraordinary ray lies again in the 
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plane of the normal and the optic axis while the ordinary ray is perpendicu- 
lar to it. 

IV. Derivation of the Electro-optic and Photoelastic Effects 

In a previous paper and in the book "Piezoelectric Crystals and Their 
Application to Ultrasonics", D. Van Xostrand, 1950, it was shown that the 
electro-optic and photoelastic effects can be expressed as third derivatives 
of one of the thermodynamic potentials. Probably the most fundamental 
way of developing these properties is to express them in terms of the strains, 
electric displacements and the entropy. For viscoelastic substances it has 
been shown that the photoelastic effects are directly related to the strains. 
In terms of the electric displacements, the electro-optic constants do not 
vary much with temperature whereas, if they are expressed in terms of the 
fields, the constants of a ferroelectric type of crystal such as KDP increase 
many fold near the Curie temperature. The entropy is chosen as the funda- 
mental heat variable, since most measurements are carried out so rapidly 
that the entropy does not vary. 

The thermodynamic potential which has the strains, electric displace- 
ments and entropy as the independent variables is the internal energy U, 
given by 

dU - Ta dSu + E m ^ + da (41) 

where Sij are the strains, 7\ ; the stresses, E m the fields, D,„ the electric dis- 
placements, @ the temperature and a the entropy. In this equation the 
strains Sij are defined in the tensor form 

Sv-tfp + p) (42) 

where the u's are the displacements along the three axis. In the case of a 
shearing strain occurring when / ^ j, the strain is only half that usually 
used in engineering practice. In order to avoid writing the factor 1/4tt, we 
use the variable 8 m = Z) m /4ir. Then, from (41), 

T - dU F - dU R - dU (AX\ 

Since, for most conditions of interest, adiabatic conditions prevail, we can 
set da equal to zero and can develop the dependent variables, the fields and 

1 "First and Second Order Equations for Piezoelectric Crystals Expressed in Tensor 
Form," VV. P. Mason, B.S.T.J., Vol. 26, pp. 80-138, Jan., 1947. 
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the stresses in terms of the independent variables, the strains and the elec- 
tric displacements. Up to the second derivatives, these are 

dE m c . 6E m 
do u dd n 



i r d 2 E 



E m _ c . 2d"-E m c . 3" -Em „ - 

as^ *'' 5 " + as^ Sij8n + aXa^ Mo . 



+ 



(44) 






2! \_i 



d 2 T,d c « , 23 2 rW ~ , d 

2! Ideas' iJqr + dS i3 d8n ijdn * db 



ndd J 



For the electro-optic and photoelastic cases, the two tensors of interest are 



d 2 Tkt 

dh n d8 

d 2 E m 



d 6 U 



= d ' En 
dSktd8 n d6o dShtd&o 



= AirmicCr, 



d 3 U 



(45) 



= (4r)f, 



95„ d5 d5 m d5 n 35 c 
For the first partial derivatives, we have the values 



dSij 



~ Gijki 



me 

dE m 



d~U _ dEn 

dS k cd8n ~ dS k C 



= 4tt# 



= —Kid 



(46) 



where C?,*£ are the elastic stiffnesses measured at constant electric displace- 
ment, hnkt are the piezoelectric constants that relate the open circuit voltages 
to the strains, and /3«„ are the impermeability constants measured for con- 
stant strain. 

With these substitutions and neglecting the other second partial deriva- 
tives, we have, from (44), 



Em = -hmijSij + D n 

Tut = fifitSu + D 



j3mn + mijmnSij + -^— D» 



h kC , nikC 
"4tt " 



wa]. 



(47) 



This equation shows that there is a relation between the change in the im- 
permeability constant due to stress in the first equation, and the electro- 
strictive constant in the second equation through the tensor m ijmn . These 
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effects, however, have to be measured at the same frequency before equality 
exists. 

To obtain the changes in the optical properties caused by the strain and 
the electric displacement we have to determine the fields and displacements 
occurring at the high frequencies of optics. Even for piezoelectric vibrations 
occurring at as high frequencies as they can be driven by the piezoelectric 
effect, these frequencies are small compared to the optic frequencies / and 
can be considered to be static displacements or strains. Hence, writing 

E m = El + £»«*", D n = D° n + D n e*", 

D = D° + D e iut , S i} = Sh 
where co = 2wf, the first of equation (47) can be written in the form 



A = -thniiSil + D° n 



/-)S r o i ?mno t-vO 

Pmn + niijmnSij + — - D c 



E„3ui r\ Jut la 8 Cl ' mno i-«0 I . ' mno r-,0 »-. ju 

m e - D n e /3 mn + »Hj mn Sij + - - - D + — D n D e 



s -| 

i*>i T\ Jul \ qS _j_ „. c I r '>™o r^O , ?mno r-vO r, jut 



(48) 



If we develop one of the fields, say Ei , this can be written in the form 

Eie*" = [0n + m.juSij + rmDl + r ll2 Dl + r llz D§D 1 e* tl 

+ [/3i 2 + ManSij + ruiDl + r li2 D\ + r l23 D° 3 }D 2 e J * 1 (49) 

+ [/3n + MijnSij + rmDl + r u2 D° 2 + r m D° 3 ]D 3 e iwi 

where the first number of r refers to the field, the second to the optical value 
of D and the third to the static value of D. Hence, for the general case, 

E m e j " 1 = D n e j "' [0 mn + m ijmn S {j + r mno D ] . (50) 

From the definition of the two tensors m i} - no and r mno given by equation 
(45), we can show that there are relations between the various components 
of the tensors. For the first tensor »i ijno , since Stj = Sji is a symmetrical 
tensor, then 

mono = Wy, no (51) 

From the definition of the tensor w,j no in the form 

4 """"° = i G-Bk) (45) 

it is obvious that we can interchange the order of 5„ and 5 so that 

"• ijno W ijon 



174 BELL SYSTEM TECHNICAL JOURNAL 

Since ij and no are reversible, it has been customary to abbreviate the tensor 
by writing one number in place of the two in the following form: 

11= 1;22 = 2; 33 = 3; 12 = 21 = 6; 13 = 31 = 5; 23 = 32 = 4 (52) 
Since the reduced tensor is associated with the engineering strains, it is 
necessary to investigate the numerical relationships between the four in- 
dex symbols and the two index symbols. From equation (48), when m 
9^ n, the change in the impermeability constant & mn is given by 

mij mn Sij-\- mji mn Sji = m T ,S r (53) 

Since S T = 2Sq = 2S# we have the relation that 

Miimn = nirsdj, m, n = 1 to 3, r, s, = 1 to 6) (54) 

In equation (45) we cannot in general interchange the order of ij and no 
since U does not contain product terms of strains and electric displace- 
ments and hence in general 

m rt 7* w 8 r. (55) 

Hence in the most general case there are 36 photoelastic constants. Crystal 
symmetries cut down the number of constants as shown in a later section. 
The tensor r mtlo defined in equation (45) as 

<*>''-"• = sSs (55) 

shows that we can interchange the order of m and n since U contains product 
terms of 8 m and 8 n . Hence 

fmno = fnmo *> ' 

and this is usually replaced by the two index symbols 

r, 10 = r mno (m, n, o = 1 to 3; q = 1 to 6). 
The so called "true" electro-optic constants are measured at constant 
strain and for this case the modifications in the impermeability constants 
are given by the equation 

E m = Dnl&n + rinoD ]. (58) 

Since m and n are interchangeable, the third rank tensor is usually replaced 
by the two index symbols 

rLo = r qo (m, n, o - 1 to 3; q - 1 to 6). (59) 

As discussed in the next sections, these constants can be determined by 
applying an electric field of a frequency high enough so that the principal 
resonances and their harmonics cannot be excited by the applied field, and 
measuring the resulting birefringence along definite directions in the crystal. 
On the other hand if we apply a static field to the crystal, an additional effect 
occurs because the crystal is strained by the piezoelectric effect and this 
causes a photoelastic effect in addition to the "true" electro-optic effect. A 
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better designation for these effects is the electro-optic effect at constant 
strain and stress. 

This latter effect can be calculated from equation (47) by setting the 
stresses Tkt equal to zero and eliminating the Sij strains. After neglecting 
second order corrections, 

E m = D n e*° l \f/L + (fine + *«*=±*S^ D °l (60) 

Since hokf/djkl = gou, the other piezoelectric constant relating the open 
circuit voltage to the stress, the electro-optic effect at constant stress can be 
written in the form 

rLo = rL„ + ™*pl!ii. (61) 

In terms of the two index symbols 

rl = rl + ^1S (62) 

4tt 

since it has been shown that g oij = g op /2 when i 9^ j, and the tensor in (61) 
has ij as common symbols which involves the summations of two terms. 

The electro-optic effect is usually measured in terms of an applied field. 
The change in the impermeability constant 0mn for this case can be de- 
termined from the first equations (47), setting Tkt equal to zero and neglect- 
ing second order terms. Multiplying through by the tensor K T op of the di- 
electric constants 

D\ = E° K T op (63) 

since the product Kl„(5l p = 1. Introducing this equation into (58) we have 
E m = ZU/C + r s mnp Kl p E° ] = D n VL + z a mno E° ]. (64) 

where the new tensor z mn0 is equal to 

~mno = 'mnp"-op • \y^J 

111 terms of the two index symbols 

4o = rf lp Ko P . (66) 

in which the repeated index indicates a summation. The difference between 
the electro-optical constant at constant stress expressed in terms of the field 
and the electro-optical constant at constant strain is 

» goij I' T S , . /S1\ 

K op = Z mno -f- niijmndpij {O/) 



4 



TV 



since the piezoelectric constants d P ij are related to the g constants by the 
equation 



rrT 



d pij = *£**£?. ( 68 ) 
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In terms of two index symbols 

£ = & + *»«<•(*. g = 1 to 6; = 1 to 3) (69) 

where a repeated index means a summation with respect to this index. 

Finally the photoelastic effect is sometimes expressed in terms of the 
stresses rather than the strains. As can be seen from equation (47), the new 
set of constants is 

ITpq = niprSrq (70) 

where the sf Q are the elastic compliances measured at constant electric dis- 
placement. 

V. Birefringence Along Any Direction In the Crystal and 

Determination of the Electro-optic and 

Photoelastic Constants 

If we take axes along the Fresnel ellipsoid when no stress or field is ap- 
plied to the crystal, the result of the electro-optic and photoelastic effects 
is to change the impermeability constants by the values 

flu = ft + Ai ; fe = ft + A 2 ; fts = ft + A 3 

#> 3 = A 4 ; fts = A 5 ; fta = A 6 



(71) 



(72) 



where 

Ai = ZnEi + Z12.E2 + ZnE 3 + wuSi -f W12S2 + W13S3 + W14S4 

+ Wi&Sb + miaSt 

A 2 = 221-El + Z22-E2 + 223^3 + TthlSl + W22^2 + W23S3 + ^24^4 

+ ithfSt + mvoSs 
A 3 = z zl Ei + Z32E2 + ZnE 3 + W31S1 + ntmS* + W33S3 + W34S4 

+ ?«35^5 + ttimSe 

A 4 = 241-El + Z42-E2 + 243^3 + MuSl + W42.S2 + tflisSs + W44S4 

+ WJ 4 frS , 6 + W4frS"6 

A 5 = 2 6 i£i + 252-E2 + 253.E3 + mnSi + W62S2 + nibsSs + w B 45 4 

A 6 = ZsiEi + Z62-E2 + 263^3 + tnnSi + MaSi + W63S3 + m*iSi 

If we transmit light along the 2' axis which, as shown by Fig. 2, makes an 
angle of degrees with the 2 axis in a plane making an angle <p with the xz 
plane, the birefringence can be calculated as follows: Keeping z' fixed and 
rotating the other two axes about z' by varying the angle $, one light vector 
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will occur when j8n is a maximum and the other when fa is a minimum. 
Using the transformation equations (31) and the direction cosines of (27), 
we find that fa is given by the equations 



,' a T 2 a -i sin 2$> sin 2^ cos . . 2 • 2 .1 
hi = Pa cos cos <p cos \f/ — (- sin <p sin \j/ 

+ 0i 2 [sin 2^ cos 2\f/ — sin 2 sin 2<p cos 2 1/' + cos sin 2^ cos 2<p] 
+ /3i 3 [ — sin 20 cos <p cos 2 ^ + sin tp sin sin 2\p] (73) 

cos sin 2<p sin 2^ 



+ j8 23 



cos' sin' <p cos' \f/ ^ "'^ °" J ** + cos 2 p sin 2 ^ 

+ /3 23 [ — sin 20 sin <p cos"' ^ — sin cos p sin 2\{/] + /3 33 sin 2 cos 2 ^ 

Differentiating with respect to \J/ and setting — = 0, we find an ex- 

pression for tan 2\J/ in the form 

— /Sn sin 2^ cos + 2/3io cos cos 2<p 
tan 2^ = + 2 fti3 sin V sin + fr»2 cos sin 2y? - 2gga sin cos g , > 
/3n[cos J cos 2 <p — sin" <p] + j8ia[(l + cos 2 0) sin 2<p] 
— /3 13 sin" ccs (p + /3 2 2(cos~ sin 2 <p — cos 2 <p) 
— /3 2 n sin 20 sin cp + /3 33 sin 2 

Inserting this value back in equation (73) we find that the two extreme values 
of fa are given by the equation 

2/3'n = 2 fa + (fcn - ft 2 )(cos 2 cos 2 ^ + sin 2 <p) + (S^ - fa) sin 2 

— /3i 2 sin 2 sin 2<? — fa sin 20 cos<p — fa sin 20 sin ^ 



± 



ten - ft 2 ) 2 (cos 2 cos 2 <p + sin 2 ^) 2 + 2(fa - fa) (fa - fa) sin 2 0X 

(cos 2 cos 2 <p - sin 2 p) + teas - #>s) 2 sin 4 - 2 ten - /3 22 )X 

[/3i 2 (sin 2<p sin 2 0(cos 2 cos 2 <p + sin 2 <p) + /3i 3 sin 20 cos <pX 

(cos 2 cos 2 <p + sinV) — & 3 sin 20 sin <p(l -f cos 2 <p sin 2 0)] 

+ 2te 33 - At) sin 2 0[/3j 2 sin 2^(1 + cos 2 0) - fa sin 20 cos <p 

4 (75) 

— fa sin 20 sin <p\ + (2/3i 2 ) 2 [sin sin 2 <p cos 2 <p + cos 2 0] 

— 4/3i 2 /3i3 sin 2 sin <p[cos 2 cos 2 93 + sin 2 <p] — 4(fafai) 
[sin 20 cos p(sin 2 ip cos 2 0+ cos 2 ip)] + (2/3i 3 ) 2 sin 2 0X 

(cos 2 cos 2 <p + sin 2 <p) — 40 13 /3 23 sin 2<p sin 4 
T + (2 fa) 2 sin 2 0(cos 2 sin 2 <p + cos 2 ?) 
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The birefringence in any direction can be calculated from equation (75) ; 
since 0n = v\/V 2 , it equals 1/ju? where m is the index of refraction corre- 
sponding to a light wave with its electric displacement in the fi' n direction. 
Similarly, for the second solution at right angle to the first, 

A -A-h ( ? « 

V y.2 
Hence if we designate the expression under the radical by K 2 and half the 
expression on the right outside the radical by Ki , we have 

i + l = A ' i; t-T-VE. (77) 

Mi M2 Mi M2 

Since mi and M2 are very nearly equal even in the most birefringent crystal, 
we have nearly 



3 



M 2 -Mi = £ = ^V* 2 . (78) 

For special directions in the crystal, the expression for K 2 simplifies very 
considerably. Along the x, y and z axes, the values are 

:j 

X, {<p = 0°, d = 90°); B x = | V(0s» - ft-2) 2 + (2ft 3 ) 2 
r, (^ - 90°, 6» = 90°); B y = | V(fci - ft*) 2 + (2/J„)» (79) 

3 

Z, (^ = 0°, 6 = 0°); B, = | V(0n - fe) 2 + (2ft 2 ) 2 - 

If any natural birefringence exists along these axes, (2ft 3 ) 2 will be very 
small compared to this and 

a 3 M 3 /l 1 \ 

b x = - (ft - fr + a 3 - a 2 ) -2^-3 + * - a.; 

B„ = f (ft - ft + A, - A.,) = f Q| - J + A, - A,) (80) 
u 3 M 3 / 1 1 \ 

Hence, for this case, measurements along the three axes will tell the differ- 
ence between the three effects Ai , A 2 and A ; , . To get absolute values requires 
a direct measurement of the index of refraction along one of the axes and 
its change with fields or stresses. This is a considerably more difficult meas- 
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urement than a birefringence measurement and requires the use of an ac- 
curate interferometer. 

If, however, the Z axis is an optic axis as it is in ADP, for example, and 
Ai = A2 = 0, a birefringence occurs due to the term fro . As shown in the 
next section, the electro-optic constants for ADP (tetragonal 4:2m) are 241 
and 2 6 3 . #63 occurs in the expression for /3 12 = A 6 , as can be seen from equa- 
tions (72), and hence the birefringence along the Z axis is 

B z = ^*2/3 12 = nlz m E,. (81) 

The constants z§z and zn have been measured independently by W. L. Bond, 
Robert O'B. Carpenter, and Hans Jaffe. Probably the most accurate meas- 
urements, and the only one published, are those of Carpenter, 5 who finds 
that the indices of refraction and the 2 6 s and z 4 i constants for ADP and 
KDP are in cgs units 

ADP 1.5254 1.4798 2.54 ± 0.05 6.25 ± 0.1 
KDP 1.5100 1.4684 3.15 ± 0.07 2.58 ± 0.05 

An even larger constant has been found for heavy hydrogen KDP by Zwicker 
and Scherrer. 6 They find at 20°C that r 63 = 6 X 10~ 7 . Using this constant, a 
half wave retardation for a X = 5461 A° mercury line occurs for a voltage 
of 4000 volts. 

For tetragonal crystals of these types the only photoeleastic constant for 
the z axis is Wm , and the birefringence for this case is given by 

B z = plmeeSt (82) 

When a natural birefringence exists for the crystal, measurements of the 
other three effects A 4 , A 5 and A 6 can be made by determining the bire- 
fringence along other directions than the Fresnel ellipsoid axes. In a direction 
of Z' lying in the XZ plane <p = 0,6 = variable and 

M 3 . /[(/3n - £22) cos 2 6 + G833 - jfe) sin 2 6 - flu sin 2 0] 2 fR ,v 

B* = 2 y + [2/S u cos 6 + 2fe sin d]\ Kfi6) 

When a natural birefringence exists, this reduces to 

cos" 6 



11 - M 



(\ - \ + At - Ao) 

+ (—-— + A 3 - A 2 ) sin 2 6 - A 5 sin 2d 



(84) 



6 "The Electro-optic Effect in Uniaxial Crystals of the Type XH 2 P0 4 ," Robert O'B. 
Carpenter, Jour. Opt. Soc. Am., in course of publication. 

6 Zwicker and Scherrer, Helv. Phys. Ada., 17, 346 (1944). 
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and hence, by measuring at 45° between the two axes, one can evaluate the 
A 5 term. 

Similarly, for the YZ plane, <p — 90°, 8 = variable and 



R _ M . A- (fin ~ fa) + (fa ~ fa) sin 2 - fa sin 2d] 2 , a _v 

n V — ~^ \ _l_ T9fl__ ™c a _ 9«._ clr, fll2 <A>/ 



•" 2 y + [20i2 cos 6 - 2/3,3 sin 0] 

Hence, when a natural birefringence exists, we have 

B yz = %\-(-2-\ + Ai - A 2 ) 

+ ( -5 - -2 + A 3 - A 2 ) sin 2 6 - A 4 sin 25 . 
\Mc fib / J 

In the XF plane 0= 90°, <p — variable and 



(86) 



2 V + [2/5i3 sin <p — fa cos <p] 2 



R _ M . /[(jSii - 0i 2 ) sin 2 *> - (fa - fa) - fa sin 2^] 2 , ? . 



(88) 



Then, for natural birefringence, 

B xy = % \(\ - A + Ax - A 2 ) sin 2 ? 
2 |_\M« Mb / 

— ( ~2 — ~z + A-, — A 2 I — A 6 sin 2<p \. 

\Mc Mb / J 

Hence, with measurements at 45° between the axes and with suitably ap- 
plied fields and strains, the three effects A 4 , A 5 and A6 can be measured. 
Since the axes of the test specimen are turned with respect to the X, F and 
Z axes, suitable transformations of the effects Ai to A6 with respect to the 
new axes will have to be made. These can be done as shown in reference (4) 
by means of tensor transformation formulae. 

Another method for measuring the constants in A 4 , A5 , A6 is to measure 
the amount they rotate the axes of the Fresnel ellipsoid. As an example con- 
sider the 241 constant of ADP. For example, if we look along theX axis and 
apply a field in the same direction, then, in equation (74), 8 = 90°, <p = and 

, r,, — 2023 —2znEi — 2^6 He 241 Ei , ,. 

tan 2^ = — = = — — — v . (89) 

fa — fa _1 2- {&> * /^AM6 — W 
2 — 2 
He Mb 

According to Carpenter, the Za electro-optic constant of ADP is 6.25 X 10 -7 
in cgs units. n a = Hb = 1.5254; n c = 1.4798; hence the angle of rotation for 
a field of 30,000 volts per centimeter = 100 stat volts cm is 

if/ = -2.25 X 10- 3 radians = 7.7 minutes of arc. (90) 
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VI. Electro-optic and Photoelastic Tensors for Various 
Crystal Classes 

Since r mnu = r nmo and z mno = z, im0 are third rank tensors similar to the 
hmij piezoelectric tensor, they will have the same components for the various 
crystal classes. For the twenty crystal classes that show the electro-optic 
effect these tensors are given below. They are given with the crystal system 
they belong to, and the symmetry is designated by the Hermann-Mauguin 
symbol. The last number of the subscript of z designates the direction of the 
applied static field. 

(91) 
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2.11 



The r tensor has similar terms. 

The photoelastic constants are similar to the elastic constant tensors 
except that m rs 9^ m sr in general. However, for the tetragonal, trigonal, 
hexagonal and cubic systems, Pockels found that Wis = IMsi . This follows 
from the transformation equations about the Z axis which is the n fold 
axes for these groups. For a rotation of an angle 8 about Z, the direction 
cosines are 



lx — — - = cos e 

dxi 


dxi 
mi = ^— = sin 

OXz 





dxi 
wi = — - = 

dx 3 


*-*5- -sin* 

dxi 


dX2 

tn* = — — cos 


e 


dx 2 

n 2 = — = 

0x3 


u = p = 

OXi 


dx 3 ft 
m 3 = — = 
dxt 




dx 3 

n 3 = r— = 1 
d^ 3 



(92) 
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Transforming the two terms mvm = mn and W2211 = m-n by the tensor 
transformation equation 



(93) 



(04) 



_ dXi dXj dXk dx( 
fflijk( — r t r r Mmnop 

ax m ox„ oXo oXp 
we find, for these two coefficients, 
m \2 — ( m n + Wit — 4w6e) sin 2 8 cos 2 6 -+- 2(w 6 2 — «u) 

sin 6 cos' 6 -+- 2(w; 6 i — w^) sin cos 6 -\- mn cos 6 -\- m»\ sin 
<m%\ = (wn + W22 — 4we 6 ) sin 2 cos 2 + 2(w ]ri — W62) 

sin 3 cos -f- 20«26 — mei) sin cos 3 6 + W21 cos 6 -\- mn sin 
If m X i = W21 for all angles of rotation we must have 

Wis + W26 = ?«61 + *W 6 2 

For all the classes that mn = m»i, either w^ 6 = — w*i6 and w; 6 2 = — ?«6i or 
else Wi6 = nu>6 = mn = m$> = 0. 

Now, if Z is a four-fold axis, as it is in the tetragonal and cubic systems, 
then, for a 90° rotation, the value of *» M or w 2 i must repeat. From the first 
of (92) this means that 

mn = W21 and W21 = mn 

For a trigonal or hexagonal system additional relations are obtained between 
mm and mn , m»2 and w E in the usual manner. Hence the photoelastic matrices 
become, for the various crystal classes, 

(95) 
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The it ten- 
sor is anal- 
ogous ex- 
cept that 
1T66 = 

7Tn — 1TI2 



From measurement 7 on the photoelastic effects at high pressure for cubic 
crystals, it has become apparent that the second derivatives of equation 
(44) are not sufficient to represent the experimental results and derivatives 
up to the fourth power should be included. This extension, however, is not 
considered in the present paper. 

VII. Photoelasticity in Isotropic Media 

The photoelastic effect in isotropic solids has been used extensively in 
studying the stresses existing in machine parts and other pieces. For this 
purpose a plastic model cut in the shape of the original is used and is loaded 
in a similar manner to that of the machine part to be studied. Since stresses 
are applied, the in photoelastic constants are most useful. If we look along 



7 H. B. Maris, Jour. Optical Society of Amer., Vol. 15, pp. 194-200, 1927. 
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the Z axis, the last of equations (79) shows that the birefringence is equal 
to 

3 

B„ = ^ V(/3i + A, - ft - A 2 ) 2 + 4(A 6 ) 2 (96) 

Since, for an isotropic substance ft = ft , we have, after substituting the 
value of Ai and A 2 , with the appropriate photoelastic constants from equa- 
tion (95), (last tensor): 

3 

B t = £ Un - tt 12 ) V(r, - T 2 y + 4zv (97) 

If we transform to axes rotated by an angle 6 about Z, the values of T n 
and T22 are given by 

Tn = cos 2 07\ + 2 sin cos 0r 6 + sin 2 6T 2 



Tm = sin 2 07\ - 2 sin cos 6T S -f cos 2 dT 2 
If, now, we choose the angle so that Tn is a maximum, we find 

+2r 6 



(98) 



(100) 



tan 20 = -J L. (99) 

1 1 — i 2 

Inserting this value of tan 20 in (98) we find 

T[ = ri ^" ^ + W{T X - Ttf + 47V 

r , = r t + r 2 _ iv / (ri _ ra)2 + 42v 

and, hence, 

7-; - r 2 ' = \/(2\ - r 2 ) 2 + 47V (101) 

Hence the birefringence obtained in stressing a material is proportional to 
the difference in the principal stresses. By observing the isoclinic lines of a 
photoelastic picture, methods 8 are available for determining the stresses 
in a model. A photograph 9 of a stressed disk is shown by Fig. 3. The high 
concentration of lines near the surface shows that the shearing stress is 
very high at these points. By counting the number of lines from the edge 
and knowing the stress optical constant, the stress can be calculated at any 
point. 

If we apply a single stress 7\ , the birefringence is given by the equation 

3 
B, = ^ (ttu - v 12 )Ti (102) 

8 See Photoelasticity, Cokcr and Filon, Cambridge University Press, 1931. 

9 This photograph was taken by T. F. Osmer. 



188 



BELL SYSTEM TECHNICAL JOURNAL 



Instead of using the constants ir u and m it is customary to use a single 
constant C given by 

B = fi e - no = r = CT (103) 

where the constant C is called the relative stress optical constant and r the 
retardation. The dimensions of C are the reciprocal of a stress and are 







■ 



Fig. 3 — Photoelastic picture of a disk in compression. 



measured in cm 2 per dyne. A convenient unit for most purposes is one of 
10 -13 cm 2 /dyne; if this is used, the stress optical coefficients of most glasses 
are from 1 to 10 and most plastics are from 10 to 100. This unit so defined 
has been called the "Brewster". In terms of the Brewster, the retardation is 



r = CTd 



(104) 



If C is measured in Brewsters, d in millimeters and T in bars (10 6 dynes/ 
cm 2 ) then r, as given by the formula, is expressed in angstrom units. 



